Radio behavior of four southern non-thermal O-type stars 
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ps| ' Abstract. We have conducted high resolution continuum observations at 

1.4 and 2.4 GHz with the Austraha Telescope Compact Array, towards the four 
southern Of stars: CD-47 4551, HD 93129A, HD 124314, and HD 150136. All 
I stars have been detected at the two frequencies. HD93129A - the only 02 I 

, star catalogued so far, and in a double system -, has also been observed at 17.8 

■ and 24.5 GHz. Its radio spectrum, complemented with previous observations at 

higher frequencies, is analyzed here. The interpretation yields the estimate of 
its mass loss rate, and a non-thermal spectral index of radiation coming from 
a putative colliding wind region. The synchrotron and corresponding inverse 
Compton luminosities are derived. 
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2 ■ !• Introduction 

. Early-type stars (O to --^02) develop strong stellar winds that strengthen while 

they are leaving the main sequence towards the Wolf-Rayet phase. The mass loss 
can be as high as 10~^ Mq yr~^ at velocities of thousands of km s^^. The winds 
^ I are optically thick, radiate through the free-free mechanism, and produce an 

. excess of flux density from the radio to the infrared range. The thermal nature 

of the radiation can be c onfirmed computing a spectral index of about 0.6 - 0.8 
( Wright fc Barlow1ll97.'j ). Once this radiation regime is established, the detec- 



tion of the wind region, and consecutive measureme nt of the flux density allow a 
straightforward determination of the mass loss rate ( Lamers Sz Leitherer 111993 ) . 



However, emission characterized by a non-the rmal spectral index has been re- 
porte d repea tedly from early-type s tars (e.g. iBieeing. Abbott. &; Churchwell I 
119891 : rBenagli ZCappa. fc Koribalskill200lh . Given the fact that the majority of 
these stars are not single, this emission is probably produced at a colliding wind 
region (CWR) of a system formed by at least two early- type stars with winds. 

The particles accelerated at the CWR can be also involved in high en- 
ergy processes. Many of these stars have been detected in X-rays (e.g. iPollock I 
|l983lj__and are prop osed as count erparts of unidentified gamma-ray sources 
((Romero~ Bcnaglia. Sz Torres ll999D . Multifrequency studies -from radio to gamma 
rays- are fundamental to completely describe the emitting regions, and help in 
a comprehensive picture of the whole stellar wind phenomenon. 
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In the frame of a southern radio survey on early- type stars, using the Aus- 

traha Telescope Compact Array (ATCA) at maximum angular resol ution, we 

have detected four non-thermal ern itters at 3 and 6 cm (Benaglia, Cap oa. Sz Koribalskil 
l2nnil : Benagli a k Koribalskil 12001 ): CD-47 4551, HD 93129A, HD 124314, and 
HD 150136. In order to better define the radio spectrum, we have performed 
ATCA observations of these sources, at 1.4 and 2.4 GHz, and here the results 
are presented. 



2. The target stars 

HD 93129A. This star belongs to the Trumpler 14 cluster, in the Carina region. 
The field is very rich in early- t ype stars, gathering 5 of the 10 earliest catalogued 



ype stars, gatnermg 

stars in the Galaxy. Following iTapia et al. I (200^ we adopt a distance to Tr 14 
of 2.8 kpc. HD 93129A is the unique representative of spectral type 2 If*. Very 
recently it was discovered to be a binary with another early-type star ( Walboriil 



I2OO2I : iNelan et al. ll200i ). 55 mas away (or 154 AU at 2.8 kpc). The star shows 
extreme values of terminal velocity (3200 it 200 km s~^), effective temperature 
(52,000±1000 K), and luminosity (log(L/L0) = 6.4±0.1) (|Taresch et al. Ill997l . 

and references therein). 

HD 93129A has been detected in the radio band at 3 and 6 cm ( Benaglia Sz Koribalskil 



I2004D . with flux densities of S'scm = 2.0 ± 0.2 mJy, Seem = 4.1 ± 0.4 mJy. The 
corresponding spectral index of 1.2 it 0.3 implies the presence of non-thermal 
emission which, in principle, is produced at a colliding wind region between the 
components of the binary system. 

HD 124314. This is an 06 V(n)((f)) field star, at a spectro-photometric dis- 
tance of 1 kpc. It has been observed three times at 3 and 6 cm with ATCA: on 
February 1998, on March 2002, and on April 2002, and presented flux density 
variations (see Benaglia, Cappa, & Koribalski 2001 for the details). It is cata- 
logued as a possible SBl bv iGies I (|l987t ). due to excursions in radial velocity 
in excess of 35 km s~^ that can be indicative of binarity. New spectroscopic 
observations are needed to confirm this result. Besides, the Washington Doble 
Star Catalogue lists a visu al companion at an angular distance of 2.8 arcsec 
( Worlev fc Douglass 



HD 150136. HD 150136 belongs to the Ara OBI association, at 1.4 kpc. This 
05 Illn(f) star has been reported as a spectroscopic binary SB2 b vlArnal et al 



(1988), with an 06 star as a close companion, in a 2.7 days orbit. iMason et al. 



(1998) listed a third component of the system, at 1.6" from the close pair. The 

ATCA 3 and 6-cm former observations detected emission from t he system at a 

level of Sscm = 2.61±0.03 mJy, and ^ecm = 5.57±0.03 mJy (BenaglTa. Cappa. Sz Koribalskil 
l200lh . 

CD-47 4551. This field star is classified as an 04 Ill(f) bv lMaiz-Apellaniz et all 



j|2003). Its spectro-photometric distance was derived as 1.7 kpc (see Benaglia, Cappa, &: Koribalskil 
I2OOI). There is no information of its binarity status in the literature. It was 
detected at 3 and 6 cm with ATCA, showing S'scm = 1-77 ± 0.05 mJy, and 
^ecm = 2.98 ± 0.05 mJy 



Non-thermal Of radio spectra 



3 



3. Observations 

The observations towards the four targets were carried out with the ATCA at 
1.384, and 2.368 GHz (20 cm, and 13 cm respectively), at the array configu- 
ration 6A, in Dec 2003. The total bandwidth was 128 MHz. The cahbration 
and analysis were performed with the MIRIAD routines. The flux density scale 
was calibrated using the primary calibrator PKS B1934-638. The resolutions 
achieved are about ~ 5" at 13 cm, and ~ 7" at 20 cm. 

Additionally, HD 93129A was observed at 17.8 and 24.5 GHz in May 2004, 
with the ATCA 6C array. The total bandwidth was 128 MHz. The flux densities 
were calibrated against Mars. The angular resolutions achieved were 2.48" x 
0.45" at 17.8 GHz, and 1.83" x 0.30" at 24.5 GHz. 

In both cases, "robust" weighted images showed the best combination of 
signal-to-noise ratio and sidelobe supression. The shortest baselines were ex- 
cluded to get rid of the diffuse emission from nearby extended sources. 

4. Results 

We detected the four targets at 1.4 and 2.4 GHz, and HD 93129A also at 17.8 
and 24.5 GHz. Gaussian fittings show all of them as point sources. The observed 
flux densities at 1.4 ad 2.4 GHz are listed in Tabic 1, together with the spectral 
indices. The flux densities detected towards HD 93129A, using the 3-mm ATCA 
receiver are: S'ly.gGHz = 1-8 ± 0.2 mJy, and S'24.5GHz = 1-5 ± 0.2 mJy. The 
spectra are displayed in Figures 1 and 2, combined with the results of previous 
observations at 4.8 and 8.64 GHz. For all detections, the polarization factors 
remain below 2%. 



Table 1. ATCA flux densities at 1.4 and 2.4 GHz 



Star 


Sp. class. t 




5'20cm 


ai3-20 






[mJy] 


[mJy] 




CD-47 4551 


04 in(f) 


3.82 ±0.20 


3.80 ±0.25 


+0.01 ±0.2 


HD 93129A 


02 If* 


7.58 ±0.50 


9.38 ±0.50 


-0.40 ± 0.4 


HD 124314 


06 V(n)((f)) 


3.70 ±0.20 


2.83 ±0.25 


+0.50 ±0.2 


HD 150136 


05 Illn(f) 


3.05 ± 0.20 


2.28 ± 0.40 


+0.54 + 0.2 



f: Mafz-Apellaniz et al. 2003 



5. The HD 93129A radio spectrum 

In close stellar systems, flux variability needs to be taken into account. An 
interpretation of the present multifrequency radio data from CD-47 4551, HD 
124314 and HD 150136 can be envisaged only once the structure of the system is 
revealed. Optical observations are currently under way to investigate this issue 
(Niemela et al., in preparation). 

We note that the component separation in the HD93129A system suggests 
a period larger than 50 yr. Here we report radio observations between 1.4 and 
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Figure 1. Radio spectra of the stars CD-47 4551 and HD 93129A 
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Figure 2. Radio spectra of the stars HD 124314 and HD 150136 



25 GHz. No flux density could be measured in the MOST 843 MHz maps due 
to source contamination. 

The relatively low angular resolution of the radio data does not allow us 
to resolve the binary system. We detect the sum of the thermal emission from 
the winds of both stars as well as non-thermal emission from the colliding wind 
region (CWR). 

5.1. Spectral contributions 

The thermal emission from the winds of both components can be approximated 
with a radio flux density of Si, oc v^'^. At the CWR, strong shocks are capable 
of accelerating wind particles to relativistic energies givin g rise to synch rotron 
emission in the presence of magnetic fields (Guinzburg Sz SvrovatskiTI Il964l. 
However, part of the synchrotron photons can be absorbed by the wind thermal 
ions. The contribution of synchrotron emission to the flux density, modifled by 
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— 2.1 

thermal absorption, can be expressed as Si, oc z^^nt gfoi' The influence of 
synchrotron self-absorption (SSA) would cause the spectrum at short frequencies 
to be represented by S^, oc We will disregard the consequences of the 

Razin-Tsytovitch effect, because according to the value of the magnetic fields 
involved, it would affect the emission up to some MHz. 



5.2. Fitting results 

We fitted the above contributions to the spectra of Fig. l:right, and found the 
best fit was given by 5",, = Ai/°-6 + Bi/^NTg-^oi^-^ i ^j^^ ^ ^ 0.17 ± 0.05, 
B = 28.6 ± 6.6, a = -1.31 ± 0.18, and tq = 1.41 ± 0.37 (see Figure 3). The lack 
of data below 1.4 GHz precluded the search for a SSA signature. 

The above expression helps us to separate thermal from non-thermal emis- 
sion, and to compute a mass loss r ate from the the rmal radio fiux. At 3 cm the 
total fiux measured is 2.0ib0.2 mJy ( Benaelia Sz Ko ribalski 200||). Consequently, 



St = 0.6 mJy, and S'nt = 1-4 mJy, and the mass loss rate for HD 93129A results 
in M = 3.6 X 10"^ M© yr'^ 

If the secondary star is an 03.5 V, like the nearby stars HD 93128 and HD 
93129B, the colliding wind regio n would be ^ 120 AU f rom the primary, and ~ 
34 AU from its companion, in an iEichler Sz Usovl ( 19931 ^ colliding wind scenario. 
A synchrotron luminosity of 6 x 10^^ erg s~^ is obtained up to 24.5 GHz. 
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Figure 3. Spectral fitting for HD 93 129 A 
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6. Discussion 

The system HP 93129A ha s a separation similar to that of WR 146 (210 AU, 
ISetia Gunawan et al. ll2nnnl 'l. The latter is a WR+0 system that has been re- 
solved in the radio continuum w ith MERLIN, and the e mission probably coming 



from a CWR could be mapped (jDougherty. Williams. &: Pollaco 2000. ). as a 38 



mas source (or ~ 50 AU at a distance of 1.25 kpc). It seems reasonable to 
assume that the size of a colliding wind region of HD 93129A should be alike, 
and thus we adopt for it a value of 40 AU. The equi partition magnetic field at 
the CWR will result in ^ 20 mGauss, (|Milev Ill98nl ). and the stellar magnetic 



field, ~ 500 Gauss. 

The same electrons involved in synchrotron processes can be scattered by 
the inverse Compton mechanism, producing gamma-ray continuum emission. 
The maximum Lorentz fac t or att ained by the electrons at the CWR, computed 
as m iBenaglia fc Romero I (I2003D . is 7max — 1-8 X 10^. This corresponds to a 
maximum synchrotron frequency of = 2 x 10^^ Hz, and a cutoff energy for 
the IC photons of E = 500 GeV. The energy at which the electron distribution 
changes, due to synchrotron and inverse Compton losses, is defined by 7b ~ 
20,000. Finally, the inverse Compton luminosity results in Lie = 1 x 10^^ erg 
s^^, well below the EGRET threshold at the location of this particular star. 

Optical spectroscopy to define the spectral type of the HD 93129A compan- 
ion, and very high resolution radio observations to resolve the system, are crucial 
to refine the figures given here. The gamma-ray predictions can be confronted 
against GLAST observations in the near future. 

In the case of CD-47 4551, HD 124314, and HD 150136, the knowledge of 
the system structure is fundamental to allow the study of the radio data. 



7. Conclusions 

Non-thermal contribution to the radio emission has been detected for the stars 
CD-47 4551, HD 93129, HD 124314, and HD 150136. The spectral fitting for the 
radio emission from the 02 If* HD 93129A enabled to extract the thermal flux, 
and determine a value for the mass loss rate of this peculiar system, of 3.6 x 10~^ 
Mq yr-^ The non-thermal radiation, coming probably from a colliding wind 
region between the stellar components, is characterized by a flux density spectral 
index of —1.3 it 0.2, and a magnetic field of ~ 20 mGauss. 
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